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Выявлены недостатки непрерывного вейвлет-преобразования не позволяющие производить их автоматизированную обработку 
при решении задачи классификации частотных составляющих сигнала. Предложен эвристический алгоритм, повышающий ин- 
формативность карт проекций изолиний вейвлет-преобразования. Показана возможность использования алгоритма для подго- 
товки информации к дальнейшей автоматизированной обработке. 


Задача исследования сигналов, в том числе, 
сигналов звукового диапазона - одна из самых по- 
пулярных прикладных задач. Часто требуется ана- 
лиз сигнала в различных представлениях: ампли- 
тудно-временном, амплитудно-частотном или ча- 
стотно-временном. Большое количество програм- 
мно и аппаратно реализованных математических 
алгоритмов позволяет адекватно переходить от од- 
ной формы представления сигнала к другой. Зача- 
стую удобные для визуального анализа формы 
представления сигнала крайне сложны для автома- 
тизированной обработки и анализа. 

Так, относительно новый математический ап- 
парат непрерывного вейвлет-преобразования 
(НВП) обладает хорошей графической интерпрета- 
цией результатов, позволяя перейти от амплитуд- 
но-временного представления сигнала к амплитуд- 
но-частотно-временному. 


Основой НВП является операция свертки сигнала 
/(?) с вейвлетом из семейства вейвлет-функций 

в ида ѵѵ, г (/) = Семейство вейвлет-функ- 

ций получают на основании материнского вей- 
влета м>(() масштабированием его параметром у и 
сдвигом на величину г: Варьируя параметрами мас- 
штабирования у и смещения г, достигают необходи- 
мой локализации как в частотной, так и во временной 
областях соответственно. Вид функции материнского 
вейвлета влияет на выявление определенных свойств 
сигнала. НВП Щт, я) определяется формулой [1]: 




Результатом НВП является двумерный массив 
чисел Щт, у) характеризующий амплитудно-ча- 
стотно-временные свойства сигнала. Значение чи- 
сла в каждом элементе массива соответствует ам- 
плитуде определенной гармоники исследуемого 
сигнала в различные моменты времени. 

На рис. 1 приведен фрагмент типичного графи- 
ческого представления НВП синусоидального сиг- 
нала в трехмерном пространстве время-частота- 
амплитуда и карты проекций изолиний. Градации 
цвета от черного к белому соответствуют измене- 
ниям амплитуды от минимальных до максималь- 
ных значений соответственно. 



Рис. 1 . Графические представления результатов непрерывно- 
го вейвлет-преобразования: слева ~ трехмерная мо- 
дель поверхности; справа - карта проекций изолиний 


Визуальный анализ результатов НВП применя- 
ется во многих прикладных и инженерных задачах 
(медицина, метеорология, оптика, и др.) [2]. Но в 
некоторых задачах исследования сигналов звуково- 
го диапазона дальнейшая автоматизированная об- 
работка результатов НВП, например, с применени- 
ем искусственных нейронных сетей (ИНС) крайне 
затруднена. Это обусловлено наличием чередова- 
ния положительных и отрицательных значений в 
массиве Щ{т, ф вдоль временной оси. Визуально 
это соответствует чередованию светлых и темных 
полос, (рис. 2) хотя, во многих случаях, темные по- 
лосы - это неинформативная часть представления. 

В данной работе предлагается ряд алгоритмов пре- 
образований результатов НВП, позволяющих устра- 
нить указанный недостаток и, таким образом, подго- 
товить их для автоматизированной обработки ИНС. 

Как известно [3], по результату операции сверт- 
ки можно судить о степени схожести двух функ- 
ций. Свойства одной из них (вейвлет-функции), 
как правило, известны, и она является базисной. 
Операция свертки является сутью НВП и опреде- 
ляет основные его свойства. 

Рассмотрим частные случаи операции свертки 
вейвлет-функций м> 8 т (1) с сигналом Д?)- 

1. Совпадение вейвлета с фрагментом функции. 

Параметры сдвига г и масштабирования 5 вей- 
влет-функции таковы, что периоды колебаний 
вейвлета Т совпадают с периодами функции Д?) 
(рис. 3 и 4). 

В такие моменты операция свертки функций да- 
ет положительное число. Причем абсолютное зна- 
чение этого числа будет тем больше, чем более точ- 
но совпадает исходный сигнал с вейвлет-функцией. 
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Рис. 2. Результат НВП сигнала вида иі)=адіп(со ] 1)+а 1 5Іп(со 2 Ѵ 



Рис. 3. Функции вейвлета ѵѵ( I) и сигнала Ц I) 



Рис. 4. Скалярное произведение Ц I) и ѵѵ( і) 


2. Совпадение инверсии вейвлета с фрагментом 
функции. 

Функции вейвлета ѵѵ(/) и сигнала /(/) в каждый 
момент времени имеют противоположные знаки, и 
скалярное произведение м>(1) и /(/) целиком лежит в 
полуплоскости отрицательных значений (рис. 5 и 6). 

В такие моменты результатом операции свертки 
функций Д/) и \ѵ(1) является отрицательное число. 
Причем, абсолютное значение этого числа будет 
тем больше, чем более точно совпадает исходный 
сигнал с инверсией вейвлет-функции. 

3. Несовпадение вейвлета с фрагментом функции. 

Функции вейвлета \ѵ(() и сигнала /(/) либо сдви- 
нуты относительно друг друга на 'А периода (в слу- 
чае совпадения периода сигнала с периодом коле- 


баний вейвлета), либо настолько различны, что ре- 
зультатом операции свертки функций/!/) и \ѵ(Г) яв- 
ляется значение близкое к нулю. 




Для решения задачи идентификации частотных 
составляющих сигнала [6], когда целью является 
выявление сигнала определенной формы и его ам- 
плитуды, значение фазы в пределах одного перио- 
да Т не играют существенной роли. Поэтому, в ка- 
честве результатов НВП можно использовать абсо- 
лютные значения НВП (рис. 7). 

В таком случае, более светлые области на гра- 
фических представлениях результатов НВП будут 
соответствовать большим значениям амплитуды 
сигнала/!/). Наиболее темные области свидетель- 
ствуют либо об отсутствии в сигнале фрагментов, 
схожих с вейвлет-функцией, либо о моментах, ког- 
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да значение сдвига г кратно четверти периода сиг- 
нала. Например, в случае НВП гармонического 
сигнала вида /(?) =а 0 $іп(<у 0 /) светлые области соот- 
ветствуют экстремумам функции, а черные - пере- 
ходам функции через ноль. 



Рис. 7 . Г оафические представления абсолютных значений НВП 

Для большинства сигналов звукового диапазо- 
на, к которым относится и музыкальная информа- 
ция [7], продолжительность звучания отдельных 
звуковых составляющих (нот) значительно больше 
периода их основной гармоники Т. Поэтому во 
временном и частотном интервалах звучания од- 
ной ноты сигнал можно рассматривать как локаль- 
но-стационарный. 

На срезе карты проекций изолиний по значе- 
нию параметра % соответствующего частоте ос- 
новной гармоники сигнала щ видно, что интерва- 
лы между пиками (всплесками), возникающими 
вдоль оси времени, соответствуют половине пе- 
риода основной гармоники этих составляющих: 
і м -іі= Т/2 и і т -ігТ, где 7М/ю 0 (рис. 8). 

В связи с тем, что наличие последовательности 
пиков на узком интервале времени говорит о нали- 
чии непрерывного локально-стационарного сигна- 
ла, то можно увеличить информативность резуль- 
татов НВП. Для этого необходимо выполнить сгла- 
живание, не потеряв информацию о наличии и ам- 


плитуде самого сигнала. Суть такого сглаживания 
заключается в формировании огибающей по эк- 
стремумам всплесков, интервал возникновения ко- 
торых соответствует рабочей частоте со и значению 
параметра 5. При изменении амплитуды сигнала во 
времени значения пиков на карте проекций изоли- 
ний будут изменяться пропорционально интенсив- 
ности всплесков каждого полупериода сигнала. 


Для формирования огибающей реализуем алго- 
ритм, который для рабочей частоты со и значения 
параметра 5 в каждый момент времени ( вычисляет 
формулу секущей Р{1) и накладывает ее на график 
колебаний амплитуд А(1) в пределах одного полупе- 
риода Т/2. 


т 


<(/ ; і 7 ' 2 ) .к/,) 

Т/ 

/2 


1 + А(^.), 


где + Т /\- 

После наложения каждой секущей, за результи- 
рующее значение сглаженного сигнала А '({) в каждый 

момент времени + /^]- берется наиболь- 

шее значение между функциями А(1) и Р(і), рис. 9. 

Далее выбирается новая рабочая частота со и со- 
ответствующее ей значение параметра ^ и операция 
вычисления и наложения секущей повторяется для 
всего временного интервала. На рис. 10 приведен 
фрагмент работы алгоритма формирования оги- 
бающей для среза на частоте со,, НВП нестационар- 
ного сигнала для трех первых семейств секущих. 

Амплитуда А{1) будет максимальна в те момен- 
ты, когда концы секущей совпадут с экстремумами 
функции А(і). То есть при совпадении частоты сле- 
дования всплесков функции А(1) (а, следовательно, 
и колебаний сигнала) с рабочей частотой со. 

В результате работы алгоритма для всех частот- 
ных срезов и для каждого момента времени полу- 
чим массив значений ]Ѵ/( т, .?) той же размерности, 
что и массив значений Щ{т, я). 




Рис. 8. Слева: карта проекций изолиний абсолютных значений НВП; справа: срез значений амплитуды по параметру 5 0 
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Рис. 9. Наложение секущих Р(і) на функцию А(і) и форми- 


рование огибающей А'( і) 



Рис. 10. Промежуточные результаты работы алгоритма для 
среза на частоте со 0 

Сформированный массив значений Щ\т, з) 
представляет собой амплитудно-частотно-времен- 


ную характеристику исходного сигнала ф{1), но, в 
отличие от результатов НВП \Ѵ[( г, з) не содержит 
отрицательных значений и чередующихся последо- 
вательностей пиков (рис. 11). 

Срез массива значений Щ\т, з) вдоль оси ча- 
стот (в любой момент времени) представляет собой 
амплитудно-частотную характеристику сигнала в 
этот момент времени. Срез вдоль оси времени - 
есть амплитудно-временная характеристика опре- 
деленной частотной составляющей сигнала (изме- 
нение амплитуды определенной гармоники на 
протяжении всего сигнала). 

Такие двумерные срезы (как и сам массив 
з)) являются исходной информацией при ав- 
томатизированной обработке ИНС, для идентифи- 
кации определенных свойств сигнала. 

Работа представленного алгоритма была опро- 
бована при обработке ряда сигналов вида 
і Д/)=а 1 8Іп(й) 1 1)+а 2 8Іп(й) 2 ?), а также фрагмента записи 
двухголосного музыкального произведения. Визу- 
альный анализ графического представления ре- 
зультатов работы алгоритма показал действитель- 
ное отсутствие в них указанных выше недостатков 
присущих результатам НВП и, следовательно, воз- 
можность дальнейшей автоматизированной обра- 
ботки. В настоящее время разрабатываются алго- 
ритмы решения задачи текущей классификации 
частотных составляющих сигнала с использовани- 
ем ИНС типа МахИеІ. 




Рис. 11. Г оафическое представление результатов преобразований значений НВП сигнала вида Ці)=а&п(і»\1:)+аіЗІп(соіі) 
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